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Abstract 

A topical collection on "Advances in Search and Rescue at Sea" has appeared in recent 
issues of Ocean Dynamics following the latest in a series of workshops on "Technologies 
for Search and Rescue and other Emergency Marine Operations" (2004, 2006, 2008 and 
2011), hosted by IFREMER in Brest, France. 

Here we give a brief overview of the history of search and rescue at sea before we 
summarize the main results of the papers that have appeared in the topical collection. 
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1 A brief history of SAR planning 



Meas uring an d pred icting the drift of search and rescue (SAR) objects has come a long way 



since 



Pingred (119441 ) made the first drift or "leeway" study of life rafts and presented it as 



"Forethoughts on Rubber Rafts". The data were unfortunately of limited value, but the ge neral 
method differed little from that of the earliest successful leeway study by IChapline who 
estimated "The drift of distressed small craft" using visual observations of drift nets to establish 
the current while simultaneously estimating the angle and speed with which the object drifted 
relative to the wind. This method of conducting leeway studies is known as the indirect method 
as it indirectly measures the moti on of the object relative to the ambient current (the leeway). 
The method reigned supreme 



Hufford and Broida 



19761 ) until the 1990s with the possible 



exception of ISuzuki and Satol ( 119771 ) who attempted to log the motion relative to the ambient 



current using a bamboo pole partly submerged and attached to the side of the ship by string. 
It should be obvious that the precision of these early experiments was not impressive, but the 
results were still of remarkable importance in the everyday work of rescue centres around the 
world. 
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In 1944 The United States Navy Hydrograph ic Office issued a manual on " Meth ods for 
locating survivors adrift at sea on rubber rafts" ( jUS Navy Hydrographic Office! . Il944l ) which 
summarized much of the current knowledge at the time of how objects on the sea surface would 
drift and how to conduct the search. The mathematical field of search theory and the wider 
topic of operations research (OR) grew out of a need to respond to the German submarine threat 
during the second world w ar. The early wo rk was pioneered by Koopman, who after having 

search and screen in g out lined the fundamentals 



provided a working manual (IKoopmanl . 119461 ) o: 
of search theory in a seminal series of papers (IKoopman 



1956a 



b 



19521 ). Without a theory of 



search the field of search and rescue would not exist and without a the ory of how the ob ject 



1980 , so 



moves, there is no way to define the search area for a moving target ( jWashburnl . 
the two fields of object drift and search theory grew up together in the post-war years. We 
refer to the combined effort of modelling the object drift and optimally allocating the search 
effort as SAR planning. In the 1950s the United States Coast Guard (USCG) first applied the 
principles of search theory to SAR planning when it published its search planning doctrine in 
a SAR manual. Since computers were not widely available, the methods were simplified and 
adapted for manual calculation. Around 1970 the USCG implemented the first computer-based 
search and rescue planning system (SARP) which was a computer implementation of the manual 
methods in the SAR manual. In 1974 the USCG implemente d the first Bayesian SAR p l annin g 
system, the Computer Assisted Search Planning (CASP), see 



Richardson and Discenza 



1980). 



CAS P was among the first applications of computer-assisted Bayesian methods (See McGravne 
for a popular a ccount of the post-war applications of Bayesian methods in search theory 
for a com p rehensive account of it s early history). For more details on 



2011 



Koopman 



1980 



and 

se arch th e ory, s ee IStond (119891 ); iFrost and Stond ( 120011 ) and the upcoming encyclopedic entry 
bv Istonel (1201.4 

CASP produced probability distributions by Monte Carlo methods, generating an ensemble 
of particle trajectories to estimate the location of the search object as a function of time. The 
trajectories accounted for the uncertainty of the initial position of the search object and moved 
the particles in accordance with a primitive drift model. This model relied on historical ship 
recordings of surface currents on a 1° x 1° monthly climatology grid and wind fields from the 
US Navy Fleet Numerical Oceanography Center (FNOC) on a 5° x 5° grid at 12 hour intervals 
forecast to 36 hours into the future. After an unsuccessful search, CASP computed the Bayesian 
posterior distribution for the location of the search object at the time of the next search by 
accounting for unsuccessful search and motion due to drift. A less coarse 3° x 3° resolution 
ocean model without tides was added in 1985. There were several evalua tions of SARP and 
CASP drift estimates using satellite tracked buoys during the early 1980s (Murphy and Allen, 



19851 ). Both SARP and CASP had mixed records at predicting the drift of search objects and 



very limited capabilities on or inside the continental shelf due to the coarse forcing fields. 

Near-real time surface current measurements near the last known position (LKP) are essen- 
tial to SAR operations. The USCG devised the self-loca ting datum marker buoy (SLDMBs) 
based on the Code-Davis drifters developed in the 1980s ( iDavisl . Il985[ ). As Argos transmitters 
became smaller and global positioning system (GPS) receivers more reliable and affordable this 
eventually led to operational use of SLDMBs in SAR operations (1 Alien! . 119961 ) . When air deploy- 
ment of SLDMBs was approved in January 2002 their use became standard routine with most 
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SAR cases, representing a major advancement in the real-time acquisition of surface currents. 
They remain an essential tool for rapidly establishing the currents near the presumed point of 
the incident. A new generation of commercially available light-weight GPS-based SLDMBs that 
can be deployed from aircraft (adhering to the NATO A-size sonobuoy standard dimensions) 
is now appearing. These new drifters have a much higher report frequency as they rely on the 
Iridium satellite network rather than ARGOS. The new generation SLDMBs will also open up 
new possibilities for physical oceanographers as the cost has come down while precision and 
reliability have improved greatly compared with earlier models. 

With the advent of high-resolution operational ocean models and the continued improve- 
ment of numerical weather prediction models (NWP), the potential for making more detailed 
predictions of the fate of drifting objects grew in the 1990s, and although the improved weather 
forecasts led to better forcing, drift models remained somewhat impervious to the advances in 
ocean modelling and numerical weather forecasting. This can perhaps best be understood in 
light of the great uncertainties in the drift properties of SAR objects. Without a proper estimate 
of the basic drift properties and their associated uncertainties, forecasting the drift and expan- 
sion of a search area remains difficult. An important change came whe n the direct method 
for measuring the leeway of a drifting object became common pr actice ( lAUen and Plourdd . 



1999 



Allen 



2005 



Breivik et al 



2011 



Hodgins and Hodginsl . Il998h . The direct method mea 



sures the object's motion relative to the ambient water using a current meter. Current meters 
small enough and flexible enough to be towed or attached directly to a SAR object started 
to become available in the 1980s, and since then almost all field experiments on SAR objects 



have employed a direct measurement technique (lAUen and Plourde 



Maisondieu et al. 



1999 



Breivik et al 



2011 



20101 ) . The direct method, together with a rigorous definition of leeway as 



Leeway is the motion of the object induced by wind (10 m reference height) and waves 
relative to the ambient current (between 0.3 and 1.0 m depth) 

and finally the decomposition of leeway coefficients in downwind and crosswind components 
makes it possible to follo w a rigorous procedure fo r conducting le e way fi eld experiments. See 



Breivik et al. 



(]201lf ) for further details. 



Allen and Plourdd fll999f ): Breivik and Alien! (120081 ): 

It was not until the 2000s that all the necessary components required for fully stochastic 
modelling using high-quality drift coefficients and detailed current and wind forecasts were 
i n place. The first opera tional leeway model to employ the USCG table of drift coefficients 
( lAllen and Plourdd . LL999J) with hig h-resolution ocean model current fields and near-surface wind 



Hackett et al. 



2006 



Breivik and Allen 



2008 



Davidson et al. 



fields went operational in 2001 (see 
2009k 

The modern era of SAR planning involving the Bayesian posterior updates after the search 
be gan in 2007 whe n USC G launched the Search And Rescue Optimal Planning System (SAROPS) 



sec 



Kratzke et al. 



( 120101 ). SAROPS employs an environmental data server that obtains wind 
and current predictions from a number of sources. It recommends search paths for multiple 
search units that maximize the increase in probability of detection from an increment of search. 
As with CASP, it computes Bayesian posterior distributions on object location accounting for 
unsuccessful search and object motion. 

By the late 2000s it was clear that although the level of sophistication and detail had grown 
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dramatically since the early days of drift nets and CASP the uncertainties in SAR predictions 
remained stubbornly high. The fundamental challenge of estimating and forecasting search 
areas in the presence of large uncertainties remains essentially the same, even though certain 
error sources have been diminished. The slow progress that has been made over the past 
decades in reducing the rate of expansion of search areas (perhaps the single best estimate of 
improvement) is an unavoidable consequence of SAR planning being at "the top of the food 
chain" in the sense tha t errors creep in from th e current fields, the wind fields, missing processes 



(e.g., wave effects, see 



Breivik and Allen 



2008 



Rohrs et al. 



£3), 



the last known position and 



not least from poor estimates of the real drift properties of the object. Indeed, sometimes 
the type of object may not even be known, effectively making the modelling exercise into an 
ensemble integration spanning a range of object categories. All these error sources accumulate 
and make SAR planning as much art as science, where rescuers still often rely as much on their 
"hunches" as on the output of sophisticated prediction tools. The fact that the majority o f 
SAR cases occur near the shoreline and in partially sheltered waters (IBreivik and Allenl . |2008) 
compounds the difficulties as the resolution of operational ocean models in many places of the 
world is still insufficient to resolve nearshore features. 



2 The state of the art of drift prediction 



Throughout the last decade these advances and obstacles to further progress have been pre- 
sented mainly through a series of workshops organized on "Technologies f or Search and Rescue 
and o ther Emergency Marine Operations" (2004, 2006, 2008 and 2011, see lBreivik and Olagnon 
20051 ) organized by the French marine research institute (IFREMER) with support from the 
Norwegian Meteorological Institute, USCG, the French-Norwegian Foundation and the Joint 
WMO-IOC Technical Commission for Oceanography and Marine Meteorology (JCOMM). As 
the last of these workshops drew near we decided that it was time to put some of the advances 
on a more academic footing by publishing a special issue, and Ocean Dynamics agreed to ar- 
range a topical collection on "Advances in search and rescue at sea". This topical collection 
focusses on recent advances in the understanding of the various processes and uncertainties 
that have a bearing on the evolution of trajectories at the sea surface, from the drift properties 
of the objects themselves to the quality of the forcing fields. 

The diffusivity of the ocean is an important factor when reconstructing the dispersion of 
particles either based on observed or modelled vector fields. In either case the dispersion is to the 
lowest order governed by the advection-diffusion equation (jTaylorl . Il92l[ ) by assuming an "eddy- 
diffusivity" coefficient. In many cases this simple stochastic m odel is sufficient for estima ting the 



De Dominicis et al 



f|2012h report 



dispersion of SAR objects over relatively short time periods, 
carefully evaluated estimates of the eddy diffusivity from a large data set of drifter trajectories 
in the Mediterranean Sea. Such regional (and possibly seasonal) estimates of diffusivity and 
the integral time scale should be carefully considered as their impact on the dispersion of SAR 
objects may be substantial. 

Stochastic ensemble trajectory models of drifting objects normally employ deterministic 
(single-model) c urrent and wind vector fie l ds and perturb the traject ories either with a random 



walk diffusivity (IBreivik and Allen 



2008 



De Dominicis et all 



20121 ) or with a more sophisti- 
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cated second-order random flight model (ISpaulding et al. 



2006 



20021 ). However, the advent of true ocean model ensembles (jBertino and Lisaeterl . 120081 ) have 



Griffa, 



1996 



BcrlofT and McWilliams 



now opened up the possi bility of exploi t ing a full vector field ensemble for estimating drift and 
dispersion in the ocean. 



Melsom et al. 



( 120121 ) compared the dispersion of passive tracers in a 
100-member ensemble of the TOPAZ ocean prediction system to the dispersion found adding 
random flight perturbations to the ensemble mean vector field and a deterministic vector field. 
The results are not conclusive in favour of the full ensemble, which is important to keep in mind 
when considering the cost-benefit of such computationally expensive operational ocean forecast 



systems. An alternative to a fu 



Rixen and 



Ferreira-Coelh 



Scott et al. 



o 



2007; 



1 model ensemble is to employ multi-mode l ensembles (see 



Rixen et al 



2008 



Vandenbulcke et al. 



2009)), which is what 



(120121 ) did when they assembled five model reanalyses and compared the weighted 



average with observ ed trajectories in the equatorial A t 



Several workers ( 



2012 



Abascal et al. 



Barrick et al 



2012 



Kohut et al 



antic. 



2012 



Frolov et al. 



2012 



Kuang et al. 



20121 ) investigated the potential for high-frequency (HF) radar monitoring 



systems to supply near real-time curre nt fields to reconstru ct the trajectories and the dispersion 



of drifting objects in the coastal zone. 



Kohut et al 



( 1201 21 ) explored the impact on search areas 



from switching to an optimal interpolation (OI) scheme for calculating total vectors from radial 



Barrick et al 



2012 



vector fields. Such techniques for extending the range of HF radars (see also 
discussed below) can make a significant difference when investigating nearshore SAR cases. 

HF radar fields and drifter studies can be used to evaluate the quality of ocean model 
current fields. Since the rate of expansion of search areas depends intimately on the quality of 
the forcing, it remains very im portant to est a blish good error estimates for each ocean model 
being used for SAR prediction. iKuang et al.l (120121 ) assessed the New York Harbor Observing 
and Prediction System (NYHOPS) using boths SLDMBs and HF currents. They found good 
agreement between model, HF radar and three drifter trajectories in the Middle Atlantic Bight 
and were able to quantify the root-mean-square differences between the modelled NYHOPS 
and the observed HF fields. 

HF short-term prediction of surface current vectors o ut to typically 12-24 hours is a tech- 



nique with great potential for near-shore SAR operations. iBarrick et al.l (120121 ) employed open 
modal analysis (OMA, see lLekien et al. 2004 ) to decompose the vector field into divergent and 
rotational modes w ithin the HF domain along the complex coastline of northern Norway (see 



Whelan et al. 



2010 



for a description of the radar deployment). They then predicted the short- 
term variation of the amplitudes of the most energetic modes based o n a relatively s hort h istory 

chose 



Frolov et al. 



of archived vector fields, giving short-term forecasts out to 24 hours, 
empirical orthogonal functions (EOFs) instead of normal modes and then employed an autore- 
gressive method to make short-term predictions out to 48 hours for an HF network in Monterey 
Bay. 

Although the direct leeway field method was established as the superior technique for estab- 
lishing the leeway of drifting objects already in th e late 19 80s, the techn i que wa s only recently 



presented in the open literature by 



Breivik et al. 



(1201 lh . 



Breivik et al. 



(j2012al ) explored how 



the technique can be applied to relatively large objects such as shipping contain ers and com- 
bined the field results with estimates from earlier work on shipping containers by 
(120021 ) to estimate how the drift varies with immersion. 



Daniel et al. 
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Most trajectory models for small surface objects ignore the direct wave excitation and 
damping since only waves whose wave l ength is comparab l e to the dimensio ns of the object will 



excert a significant force on the object ( jBreivik and Allen 



2008 



Mei, 



19891 ). Since SAR objects 



are typically smaller than 30 m their resonant ocean waves will have only negl i gible energy 
Howe ver, waves will also affect an object through the Stokes drift ( Phillips . 19771 Holt huii sen . 



2007), which is a Lagrangian effect not visible in an Eulerian frame of reference. iRohrs et al. 



(120121 ) explored how the Stokes drift affects surface drifters with and without leeway directly 
and through the addition of the Coriolis-Stokes effect to the momentum equation. The term 
adds an additional deflection to upper-ocean currents caused by the Coriolis effect acting on 
the Stokes drift. This has clear relevance for the operational forecasting of SAR objects as well 
as for the interpretation of SLDMB trajectories, although it is not clear yet how large the effect 
is for real-world search objects that also move under the direct influence of the wind. 

Finally, the importance of being able to estimate the point of an accident based on a 
debris field was made poig nantly clear after the AF447 aircraft accident on 1 June 2009 in 
the equatorial Atlantic (see 



Stone et al. 



2011 



for an account of the search effort following the 
accident). Using SAR trajectory models for backtracking is not trivial since it effectively 
means reversing the (usually weakly nonlinear) processes that propel the object. In principle 
it is better to run a model forward and iterate, as lBreivik et al.l (I2012bl ) has demonstrated, but 
nevertheless direct ba cktracking can be employed if the model integration times are modest. 
Drevillon et al.l (120121 ) describes the amount of preparation that went into the so-called "Phase 



III" of the search. Detailed regional atmospheric reanalyses and ocean model hindcasts were 
performed to prepare a multi-model high-resolution ensemble of wind and curren t fields that 
were t hen used to perform a range of backtracking trajectory integrations. Similarly, IChen et al. 



(]2012l ) included a wind drag factor and were able to estimate the point of impact for the 
AF447 accident based on backtracking the observed debris field. The method of using a w ind 
drag coefficient to fine-tune the drift properties was also employed by lAbascal et al.l (120121 ) to 
investigate the optimum balance of HF current fields and wind fields required to backtrack 
drogued and undrogued drifters. 

The 12 articles in this topical collection provide a snapshot more than a complete overview of 
the state of object drift modelling and SAR prediction at sea as it stands today. We hope that by 
putting together this special issue we provide a starting point for new workers in the field as well 
as a body of references of what has been published earlier. This is particularly important in an 
operational field such as SAR planning where a majority of the work to date is "grey literature" 
in the form of technical reports that may not be readily accessible or properly vetted through 
peer review. SAR planning and object drift modelling demand both mathematical rigour and 
experimental finesse to advance further. Peer-reviewed communication is the most efficient 
way to achieve this. It is our hope that this special issue will contribute to a more academic 
approach to this exciting field. 
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